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Abstract 

Fission fragment mass and angular distributions were measured from the 
64 Ni+ 197 Au reaction at 418 MeV and 383 MeV incident energy. A detailed 
data analysis was performed, using the one-body dissipation theory imple- 
mented in the code HICOL. The effect of the window and the wall friction 
on the experimental observables was investigated. Friction stronger than one- 
body was also considered. The mass and angular distributions were consis- 
tent with one-body dissipation. An evaporation code DIFHEAT coupled to 
HICOL was developed in order to predict reaction time scales required to de- 
scribe available data on pre-scission neutron multiplicities. The multiplicity 
data were again consistent with one-body dissipation. The cross-sections for 

touch, capture and quasi-fission were also obtained. 
25.70.Lm 
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I. INTRODUCTION 



The mechanism through which energy is dissipated in low energy heavy ion reactions is 
a topic widely discussed in the last 20 years. The study of deep-inelastic collisions in the 
1970s stimulated the development of dissipative theories of nuclear dynamics Following 
this, the quasi-fission phenomenon was discovered in the 1980s. As well as fusion-fission, it 
belongs to the class of the completely damped reactions, where all the initial kinetic energy 
is dissipated, but depending on the entrance channel, the mass asymmetry degree of freedom 
may be fully or partially equilibrated . Large mass transfer is achieved in a relatively short 
time ( 2-10 xl0~ 21 s) [|J and the system quickly re-separates due to the absence of a fission 
barrier. 

In heavy symmetric systems, quasi-fission is the dominant reaction mechanism. A num- 
ber of systematic studies done at GSI |§-|6| using 238 U and 208 Pb beams incident on targets 
ranging from 16 O to 89 Y, have identified the following experimental signatures of quasi- 
fission: 

• Fragment mass distributions wider than the mass distributions resulting from fusion- 
fission reactions. 

• Asymmetries in the mass-angle correlations increasing with the target mass. 

• Angular anisotropies significantly larger than those in fusion-fission reactions. 

Because of their unique nature, quasi-fission reactions can provide a testing ground for 
different models of dissipative heavy ion collisions. While in fusion-fission reactions, the 
fragment mass distributions are not informative of the reaction dynamics, in quasi-fission, 
the mass drift towards symmetry is a sensitive probe of the dissipative forces acting between 
the reaction partners. The asymmetry in the mass-angle correlations, wherever present, is 
directly related to the time the system spends in the coalesced state, before re-separation. 
In the GSI measurements it was determined, that for all targets heavier than 40 Ca the sys- 
tem made less than one turn before re-separating and the turning angle of the intermediate 

2 



complex was used to extract the dynamical time scale of the reaction. The data analysis 
found that the relevant friction mechanism is one-body in nature and is temperature inde- 
pendent. As a result of these works, the extra-push model parameters were established in a 
systematic way. They were also used in the development of the one-body dissipation model 
implemented in Feldmeier's code HICOL J7J. 

Another experimental approach for determining the reaction time scale is the measure- 
ment of pre-scission neutron multiplicities. This approach is not limited to quasi-fission 
reactions, but was also used to deduce fusion-fission time scales HU. When the excitation 
energy of the coalesced system is high enough, the neutron evaporation times are much 
shorter than the scission time and thus, the pre-scission neutron multiplicity can be used 
as a clock. The theoretical treatment of the data involves a statistical model code, which is 
used to predict the pre-scission neutron lifetimes and to deduce the reaction time. Although 
the dynamics of the reaction cannot be treated with this model, and a direct comparison 
with the different dissipation theories has not been done, the quasi-fission time scales de- 
rived from pre-scission neutron multiplicities are comparable to the results obtained with 
mass-angle correlations in similar systems and temperature regimes || [| . 

In a recent work, Wilczynski et.al. []rO|,|TT] re-analyzed Hinde's et.al. [RLRJ pre-scission 



neutron multiplicity data, using an evaporation code DYNSEQ coupled to the outputs of 
HICOL. In that way, the dynamics of the reaction was taken into account. The time scales 
obtained were an order of magnitude longer than the original results of Hinde et.al. Sur- 
prisingly large values for the dimensionless friction coefficient 7 = 17 — 50 were reported. 
These were derived after scaling the one-body dissipation in HICOL by a factor k s = 4 — 12. 
Using this new analysis, and results of analysis |12) of the GDR 7-rays in coincidence with 
fission for the O+Pb, S+W,Pb reactions [DJ, the authors claimed to observe strong two- 



body dissipation setting in at nuclear temperature of about 2 MeV, with 7 following 1/T 2 
temperature dependence. One-body friction, for comparison, yields 7 ~ 4 independent of 
temperature. 

The strength of two-body dissipation is still an open question in the literature. While 



most authors find that two-body dissipation is weaker than one-body and yields shorter 
pre-scission time scales fl~4] , |i~5|| , there are some theoretical works that predict strong two- 
body viscosity and long time scales associated with this (e.i. |16[). Such strong dissipation, if 
present, would influence significantly the quasi-fission reaction dynamics and the observables 
associated with it: fragment mass distributions and mass-angle correlations. It is desirable 
to compare dissipation results deduced from pre-scission neutron multiplicities and from 
fragment mass and angular distributions. At present, there are four quasi-fission reactions 
in which the pre-scission neutron multiplicities have been measured ||. The data set of mass- 
angle correlations is larger, but until now, there was no measurement of both observables 
in the same reaction at the same center of mass (cm.) energy, and hence, the same nuclear 
temperature. 

In this work the mass and angular distributions were measured in the 64 Ni+ 197 Au reac- 
tion at the 418 MeV and 383 MeV laboratory energies. The mass-angle correlations were 
also constructed. For the first bombarding energy, pre-scission neutron data exist and a 
comparison of the two experimental approaches is possible. The second energy was used to 
investigate the temperature dependence of the friction. With the development of the code 
HICOL, a detailed analysis of the mass distributions became possible and the relative con- 
tribution of the window and the wall friction could be determined from the data. An upper 
limit to the reaction time, and the friction, was determined by an analysis of the mass-angle 
correlations using the one-body dissipation model HICOL. An evaporation code DIFHEAT, 
coupled to HICOL, was developed and applied for analysis of the pre-scission neutron data 
published in Ref. and a comparison with Wilczynski's results (obtained with a different 
code) is presented. The deep inelastic, capture and touch cross-sections were also obtained. 

II. EXPERIMENT 

The experiment was performed at the Stony Brook Nuclear Structure Laboratory, using 
64 Ni beam from the FN Tandem Van de Graaff and Super-conducting LINAC accelerators. 
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The reactions, for which the pre-scission neutron multiplicities were measured and the 
two existing interpretations for the reaction times, determined in Refs || and QTUJ are given 
in Table |. To complement the pre-scission neutron multiplicity data, we measured the 
mass and angular distributions of the fragments produced from 64 Ni + 197 Au at 418 MeV 
incident energy. Although, in this case the discrepancy shown in Table | is mild compared 
to the 64 Ni+ 238 U and 40 Ar+ 238 U reactions, the 197 Au target was chosen for its lower fissility. 
The measurements were performed in "singles" mode with Silicon Surface Barrier detectors 
(SSBs) and the event reconstruction was done assuming binary fission. 

The experimental layout in the 2.4 m diameter "BigMac" scattering chamber is shown 
in Fig.|l| Eight large area SSBs (EG&G Ortec BF-028-400-60) were mounted on a movable 
platform 10° apart located at distances 50 cm and 40 cm from the target. Measurements 
were done at two positions of the platform, covering from 20° to 90° laboratory angles in 
the forward position and from 80° to 150° in the backward position. The platform was 
connected to a cooling system and the detectors were cooled to — 10°C, which allowed over- 
biasing without damaging the detectors. Operating the SSBs in high field strength was 
crucial, since this minimizes the pulse-height defect and the plasma delay. For the heavy 
ions studied in this experiment, the plasma delay time was between 0.2 — 1.5 ns. Corrections 
were made to compensate for this effect. The pulse height defects were E p hd = 5 — 25 MeV. 
Although systematic studies of the two effects listed above are present in the literature 



17| ) [18| , the correct mass reconstruction for ions as heavy as the fission fragments measured 
in this experiment was difficult. An original procedure for mass reconstruction was developed 
in this work and is discussed in Ref. fL9|| . In addition, two small area ion-implanted silicon 
detectors (EG&G Ortec BU-013-25-300) were placed at ±20° with respect to the beam axis 
at a distance 81 cm from the target to monitor the beam quality and for normalization of 
the measured fission fragment yield to the Rutherford scattering cross section. 

The reconstruction of the mass and the energy of the detected fragments, requires an 
energy and absolute Time-Of-Flight (TOF) calibration of the detectors and the subsequent 
electronics. The energy calibration was done using a 228 Th a- source and the elastic peak 



from a 247 MeV 64 Ni beam. At the Ei a b = 418 MeV energy, additional calibration points 
from elastic scattering were used for the detectors that were within the reaction grazing 
angle. The energy resolution was ~ 2% for the elastically scattered beam and ~ 1.2% for 
a-particles. An absolute timing calibration was done using the reference RF signal of the 
LINAC accelerator and the elastic peak from the 247 MeV run. A thin (100 /ig/cm 2 ) self- 
supporting Au target was used in the calibration and a 250 /ig/cm 2 Au target was chosen 
for the measurement, compromising between good timing and a reasonably high count rate. 
The width of the beam pulse at the target was < 600 ps FWHM. Corrections were made 
for the plasma delay time of the heavy fragments, and the overall timing resolution for the 
elastic reaction products measured was < 1 ns. 

In addition, the relative solid angles of the SSB's and the monitors were calibrated using 
the Rutherford scattering yield. Several calibration runs at 247 MeV, bracketing in time the 
higher energy data, were recorded during the course of the experiment. 

III. DATA ANALYSIS 

A. Mass Reconstruction 

Two-dimensional (2D) scatter plots energy versus TOF were made for each detector, 
and the reaction products were separated from the slit scattered beam and low-energy- 
noise backgrounds by applying 2D gates. An example of such a plot and the gate applied 
for detector 3, located at 40° in the lab system, is shown in Fig.^(a). The high energy 
elastic peak dominates the yield. Quasi-elastic events are seen at an energy lower than 
the elastics and with similar TOF. The regions of the deep inelastic and the quasi-fission 
fragments are indicated in the figure. Although the elastic and quasi-elastic channels are 
not of interest for the quasi-fission reaction studied here, these events (wherever present) 
provided an important check for the mass reconstruction, so they were included in the first 
step of the data analysis. Later, when only the damped reaction products were studied, 
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total kinetic energy (TKE) versus mass scatter plots were made for each detector, and a 
gate was applied to select the deep inelastic and quasi-fission events only. See Fig. 0(b) for 
6 = 40° (lab). 



After selecting the events of interest, the mass was reconstructed, as described in Ref |L9 
with overall mass resolution better than 10 mass units. 



B. Efficiency Simulations 

A Monte-Carlo simulation program was developed to determine the detection efficiency 
for each mass at each detector position. Events were generated in the cm. with flat mass 
distribution and da/dQ c . m . = 1/ sin 8 (assuming isotropic distribution in the polar angle 8). 
The kinetic energy of the fragments was drawn from a Gaussian distribution with parame- 
ters determined from the measured TKE distributions, as discussed later. The simulations 
included energy loss in the target and in the front Au electrode of the detector, pulse- 
height-defect and plasma delay, electronics energy thresholds, TOF and velocity cuts. After 
applying the experimental cuts and the experimental mass reconstruction procedure, the 
events were transformed back into the cm. and the "detected" yield was compared to the 
input. The efficiency curves for the two energies studied and for the two positions of the 
platform holding the detectors are shown in Figure || For small masses, the input angular 
distribution is recovered from the "data", in the intermediate mass region (A= 125 — 130), 
the detection efficiency is reduced in the three most backward detector positions, and the 
heavy fragments (A> 150) can only be detected in the forward position of the platform 
(from 20° to 90° in the laboratory system). 

The mass cuts observed in this experiment agree very well with the simulated efficiency 
curves. This gives confidence that the mass reconstruction in the experimental data was 
done correctly, and also, that all relevant effects were included in the simulations. Based 
on the simulated efficiency, correction to the experimental yield was done for the mass bins 
in which the efficiency was between 0.5 and 1. For lower efficiencies, the centroid of the 
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fragment kinetic energy distribution is not covered, so these data were excluded from the 
analysis. 



IV. RESULTS AND DISCUSSION 

A. Fragment Mass-TKE Distributions 

2D scatter plots of TKE (in the cm.) versus Mass were generated for each detector. The 
data were corrected for pre-scission and post-scission neutron emission, using the measured 
quantities from Ref. ||: for 418 MeV, v vre = 3.15, v post = 2 — 5 depending on the mass of 
the fragment. For E; a fe = 383 MeV, experimental data are not available, so values estimated 
on the basis of HICOL's predictions for the heat along the trajectory and in each fragment 
at the scission point were used: v pre = 1.5 and v vost = 1 — 3 increasing with the fragment 
mass. In Fig. |2|(b) an example of TKE-Mass scatter plot is shown. The line indicates the 
Viola-systematics energy. The measured centroids of the TKE distributions agree with Viola 
systematics and with previous quasi- fission studies The variance of the distribution, 

for which systematic data and empirical formulae do not exist, was determined to be m 15% 
from the centroid and is again consistent with the measurements in Ref. ||. The HICOL 
calculations, discussed later, predict correctly the centroid of the TKE distribution, which 
shows that the friction mechanism implemented in the code produces the right amount of 
damping. 



B. Fragment Angular Distributions 

1. Experimental Results 

The energy-integrated cross-sections binned in 10 mass units, = 2tt 

were constructed, after lab to cm. transformation and normalization to the Rutherford 
scattering yield. 
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Figures |4] and [5] show the results from the two energies measured. The data for the 
masses up to the symmetric mass split are shown with full points. Above symmetry, only 
a limited angular range was covered, due to the various factors restricting the measurement 
and the mass reconstruction. The data obtained in this region were used for a consistency 
check and are plotted with open points in their complementary mass bins. The error bars 
shown in the figures represent the statistical errors from the measurement. The data show a 
significant forward - backward asymmetry typical for short-time scale reactions in which the 
rotational period of the system is longer than the reaction time. This behavior is consistent 
with one-body friction. 



2. HICOL Calculations Along the Trajectory 

A series of HICOL calculations was performed to investigate the dissipation mechanism 
responsible for the forward-backward asymmetry in the experimental angular distributions. 
The code implements a macroscopic description of dissipative heavy ion collisions |7|] . A set 
of collective shape and angular coordinates is used and their evolution is followed, solving the 
equations of motion. The coupling of the macroscopic and the intrinsic degrees of freedom is 
considered explicitly in the derivation of the friction and the mass diffusion tensors, in which 
Fermi gas one-body formulae are used. A typical "along-the-trajectory" calculation, which 
will be used as an illustrative example is shown in Fig. || We will consider a quasi-fission 
reaction of 64 Ni+ 197 Au at 418 MeV laboratory energy and initial orbital angular momentum 
L= 50h. The potential energy, the kinetic energy in the shape degrees of freedom T s , in 
the angular degrees of freedom (T re i and Tj nt ) and the dissipated energy are plotted as a 
function of time. In panel (a), the actual shapes along the trajectory are also shown. The 
trajectory starts at time t — when the nuclei are 14 fm apart. All the energy available to 
the system resides in kinetic energy of the shape and relative rotation. Within 6x 10 _22 s, the 
nuclei come into contact. At that time, the proximity friction has already produced about 
8 MeV of heat. After the touching point, as the window opens, the radial motion is rapidly 
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damped. With this, the shape loses most of its kinetic energy, since the other two shape 
coordinates - the neck and the asymmetry do not possess any initial momentum. The viscous 
drag in the window sets into motion the intrinsic rotation. At t — 15xl0~ 22 s the window has 
produced 68 MeV of heat and is already widely open. At this point the shape changes from 
a necked-in to a convex type. The system enters the mono-nuclear regime, and a uniform 
wall is considered for the subsequent evolution. A long creeping motion with practically no 
kinetic energy in the shape follows. Without the window friction against the net particle 
flux, the asymmetry degree of freedom begins relaxing. It has been severely hindered up 
to now. As we move along the trajectory, the deformation develops again, the window 
comes back and the total dissipation increases. The potential is rapidly falling, the system 
gains back kinetic energy in the shape (10 MeV in this case) and it scissions separating 
with final angular momentum Lf = 43h and final fragment mass numbers Ai = 137 and 
A 2 = 124. The TKE of the fragments is 221 MeV, consistent with our measurement and 
with Viola systematics. The mass variance is a AA = 888 and the total number of exchanged 
particles through the window - N ex = 6614. Within the reaction time of ~ 140 x 10 _22 s the 
system does not complete one full rotation ( the scattering angle is —54°). The fragment 
masses are correlated with the emission angles, which will result in asymmetric mass-angle 
distributions. This is a typical quasi-fission trajectory, where the system scissions bypassing 
the compound nucleus stage. The mass asymmetry is not fully relaxed, but almost all the 
initial kinetic energy is dissipated and the fragments fly apart driven by the strong Coulomb 
repulsion. 

In the model implemented in HICOL, there are no free parameters. The only ad hoc 
procedure used is in the combination of window-plus- wall and mono- wall friction. For shapes 
with a "fat" neck, there is an ambiguity in the choice of a window plus two walls, or a mono- 
nucleus shape without a window. To make a smooth transition between the two regions, a 
form factor f(x n ) multiplying the rate of energy dissipation Q, is introduced as follows: 

Q = (1 - f(x n ))Q win+wal1 + f(x n )Q mon °, (1) 
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with x n = R nec k/min(Ri, R2) measuring the ratio between the neck radius and the radius 
of the smaller nucleus. For "thin" necks (0 < x n < 0.8), f(x n ) = and the full window- 
plus-wall friction is used. In the region 0.8 < x n < 1, where the shape makes the di-nucleus 
to mono-nucleus transition, f(x n ) is smoothly rising from to 1 and both widow-plus-wall 
and mono- wall are used. For x n > 1, the energy dissipation is calculated with the mono- wall 
formula. In the outgoing part of the trajectory, when the neck develops again, f(x n ) drops 
from 1 to 0, after which, the full window-plus- wall friction is used until the scission point. The 
functional form of f(x n ) and the limiting value x n = 0.8 are arbitrarily chosen. The concern 
when using such a description comes from the fact that the two formulae ( window-plus- wall 
or mono-wall) give different results, when applied to the same shape. The window-plus-wall 
formula can produce 4 to 8 times stronger friction against the mass- asymmetry equilibration 
than the wall formula alone. Fig. |^ (a) illustrates the mass evolution of the projectile-like 
and the target-like fragment along the = 418 MeV, L= 50 h trajectory. The evolution 
of the form factor f(x n ) is shown below it in Fig.[7| (b). One can easily see, that the main 
mass transfer from the target to the projectile happens during the mono-nucleus part of the 
trajectory, when the window friction is turned off. In Fig. [7] (c) and (d) we have calculated 
the mass evolution in the two extreme cases: mono-wall friction along the whole trajectory, 
or full window-plus-wall everywhere, respectively. Both these cases are clearly unphysical. 
Without the window friction, all trajectories proceed rather rapidly to symmetry, which is 
in contradiction with the experimental results (Section [IV D|) , while with the window-plus- 
wall acting everywhere, there are no trajectories that produce symmetric mass-splits. The 
realistic description of the reaction dynamics needs a mixture window-plus-wall and mono- 
wall friction. We consider that the form factor f(x n ) should be used as free parameter and 
adjusted to match the experimental mass distributions. 

In the dynamical analysis of Hinde's et.al. pre-scission neutron data |9|], Wilczynski, 
Siwek-Wliczynska and Wilschut |T]J scaled the energy dissipation rate Q by an additional 
factor k s , which was treated as a free parameter. The factor k s was allowed to obtain values 
smaller than unity, when a reduction of the full one-body friction was needed, and k s > 1 

11 



was interpreted as evidence for two-body dissipation. Along the trajectory, the authors of 
Ref. ]TI|] divided k s into k l s n and k° ut . In the early stage of the reaction, when the system is 
reasonably cold, a factor /c* n = 0.5 was used, and after the mono-nucleus stage is reached and 
the system has converted a significant amount of kinetic energy into heat - k° ut = 4 — 12 was 
applied to make the reaction time long enough to allow the emission of the requisite number 
of pre-scission neutrons. The justification for using such a description was the expected 
temperature dependence of nuclear viscosity. 

Fig.| shows the L=50/i trajectory calculation for 64 Ni+ 197 Au at = 418 MeV, which 



was done using k l s n = 0.5 and k° ut = 4 as described in Ref. [ITJ. The quantities plotted are 
the same as in Fig|| which was for the standard (k s = 1) HICOL trajectory. By comparing 
the two calculations, one sees that the shape of the potential energy, kinetic energy and total 
heat curves is similar, but the time scale is stretched in the large friction trajectory. The 
relative contributions from window and wall dissipation are different. By scaling down the 
friction in the early stage of the reaction when the window-plus-wall formula is used, the 
window contribution is reduced. Then, in the mono-nucleus stage, the weaker wall friction 
is scaled by a factor of 4. This brings the system into a relatively flat region of the potential 
with very little kinetic energy (mainly in rotation) and strong friction hindering its motion 
towards scission. It takes 3.5 times longer to develop the deformation and find its way to 
the steep slope of the potential that leads to scission. During this long reaction time, the 
system makes more than two full rotations, as seen in the shape evolution in the figure. 
The calculated final scattering angle is —760°. This will lead to orbiting type (~ 1/ sin#) 
angular distribution of the fragments, instead of the asymmetric distribution expected from 
the trajectory in Fig.|| which agrees qualitatively with the experimental results. The mass 
drift in the large-friction trajectory is similar to the one-body case, since the window friction 
was initially reduced (k l s n = 0.5), and later increased, in the final part of the trajectory. 
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C. Cross-Sections 



The experimental angular distributions were extrapolated to 0° and 180° as described 
by Back et.al. in Ref. [PPf] . For fusion-fission reactions the angular distributions of the 
fragments are given by [ETT : 



1=0 

2 



x- 



Here I is the total spin of the fissioning nucleus, K is the projection of / on the fission 
axis, Tf is the transmission coefficient for fusion of the 7th partial wave and d QK {9) are the 
real parts of the Wigner Z)-functions. Eq. ^| gives angular distributions symmetric about 
90° which approach the 1/ sin# classical distribution for I » 0. These are not applicable 
to describe the asymmetric angular distribution in short-time scale quasi-fission reactions. 
However, since a theoretical prescription for this type reactions is not available, Back et.al. 
20| used the following expression to fit the data from a measurement similar to the one 



done in this work: 

da 



2tt sin 6 c . m . exp((3(6 c . m . - tt/2)) W{9 c . m ), (3) 



dQc.rn. 

where j3 is the angular slope parameter. The same approach was adopted here. The trans- 
mission coefficients Tj were calculated using the extra-push estimate to the complete fusion 



cross-section with a program due to Back |22|. The Wigner D-functions for I, K < 25 were 



were 



calculated using the CERN library routine DDJMNB, and the recursion relations |23| 
used for larger spins. The value of K was determined from the mass bin closest to symmetry, 
and kept fixed in the other mass bins. Best results were obtained with = 25t\o f° r the 
418 MeV data and Kq = 15^5° at 383 MeV. The overall normalization and the parameter (3 
were varied to produce best fits to the data. In the systems studied by Back et.al. PU| the 



clean separation of the quasi-fission and deep inelastic events was possible, since the mass 
distributions were peaked at symmetry. In the case studied here, the mass distribution has 
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a minimum at symmetry and these two reaction channels could not be separated (see also 
Fig|| The deep inelastic scattering exhibits a maximum at the grazing angle of the reaction 
for A < 80 (right panels in Figures |] and [5]) and the above description of the angular 
distribution does not apply The laboratory grazing angles for Ej a & 418 MeV and 383 MeV 
are 55° and 65°, respectively. In order to fit the data, a Gaussian distribution was added to 
the distribution given in Eq.[3] and the relative contributions of the two distributions were 
adjusted to obtain agreement with the measurement. With this, an attempt was made to 
estimate the deep inelastic and the capture cross-sections separately. For the A > 80 mass 
bins (left panels in Figures |] and |5]), the cross section is only due to capture reactions. The 
values for the total cross section measured, a toU ch = &di + ^capture, and the estimated con- 
tributions from the two reaction channels are given in Figs. 4, 5. The errors in a toU ch, around 
25%, come from the quality of the fits and the uncertainty in determining the monitor po- 
sition and the detector solid angle ratios. The contribution from the statistical errors of the 
data is much smaller. The upper limit of the deep inelastic cross section <jdi was estimated 
with the assumption that for A < 80 all the cross-section forward of the grazing angle is 
in this channel. This brings the corresponding reduction in the capture cross- sect ion. The 
angle-integrated cross-sections for the different mass bins are listed in Table |T[ 



D. Mass Distribution 

After integrating the angular distributions in the different mass bins, the mass distribu- 
tions for the two energies measured were obtained (see Fig. H). At both energies, the dis- 
tributions differ dramatically from the mass distributions of fusion-fission reactions, which 
are narrower and peaked at symmetry. This behavior can be qualitatively explained since, 
for this system which has fissility x = 0.915, the symmetric saddle does not exist. The 
potential energy landscape always favors asymmetric mass division. Theoretical mass dis- 
tributions were obtained using HICOL's predictions for the final mean fragment masses and 
their variances for each trajectory that leads to a damped reaction. A weighted sum over 
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angular momentum was constructed, using the measured cross-sections. The resulting mass 
distributions are plotted together with the experimental data in Fig. |9|. The dashed-line 
histograms give the HICOL predicted mass-distributions. Although the general trend of the 
data is reproduced, a quantitative agreement is not achieved. The code under-predicts the 
deep inelastic scattering and over-predicts the yield at the symmetric mass division for both 
energies. Since the mass drift is governed by the window friction, and in HICOL the relative 
contributions of window-plus-wall and mono-wall friction are mixed with a form factor deter- 
mined from the shape of the system along the trajectory, as discussed above, a variation of 
the form factor was considered in this work, in order to reproduce the experimental results. 
Allowing for the window friction to contribute for larger R nec k by changing the limits of x n 
from 0.8 < x n < 1 to 0.95 < x n < 1 gave a satisfactory result for the 418 MeV data and a 
trend in the right direction for the 383 MeV data (the solid-line histograms in the figure). 
Further increase in the lower limit of x n still did not reproduce the large deep inelastic cross 
section measured at Ei ab = 383 MeV. 

The discrepancy between HICOL and the data in the deep inelastic channel was ad- 
dressed by Feldmeier and was attributed to the fact that the narrowly-peaked phase 
space approximation is not strictly valid for the deep inelastic trajectories. The mass vari- 
ance is under-predicted in this region, which limits the partial waves that can contribute 
to the cross-section in the near-projectile and target mass bins. Including fluctuations in 
the trajectory model will probably improve the agreement with the data. The code gives 
a reasonable description of the 418 MeV data, in which the deep inelastic events comprise 
~ 1/4 of the total touch cross-section. At the lower energy, where almost half of the touch 
cross-section is in the deep inelastic channel, only the shape of the distribution is reproduced 
by the calculations. 

Comparison of the = 418 MeV mass distribution data with other theoretical de- 



scriptions from the literature is presented in Fig. |10| In the work by Shen et.al. ||, it was 
found that the full window-plus-wall dissipation gave reaction times which were too long 
and inconsistent with the time scales derived from the mass-angle correlations. The authors 

15 



came to the conclusion that the window friction should not be included in the calculation. 

In Fig. [TO) (a), a calculation without the window friction is presented. It differs dra- 
matically from the data, since the window friction is the main mechanism that prevents 
the fast equilibration of the mass asymmetry. Fig. [10] (b), shows a calculation done using 
Wilczynski's (JTXJ] prescription for scaling the friction along the trajectory, which was already 



discussed in Section [IV B 2[ Although this calculation reproduces the mass distribution, it 



does not reproduce the mass-angle correlations, since the reaction times are longer than the 
rotational period of the system (See Section [I V B 2|) . 

A frequently used approach to implement friction in various descriptions of heavy-ion 
reactions is to use window-plus-wall dissipation during the whole reaction, but scaled by a 
factors 0.3, which was determined from fitting experimental fission fragment TKE distribu- 
tions and giant resonance widths [ 24 1 . Figure [TpJ (c) shows the calculated mass distribution 



using this approach, in comparison to the 418 MeV data. With the window friction present 
during the whole trajectory, even if it is scaled down, the mass drift towards symmetry is 
hindered significantly, resulting in disagreement with the data. 

Using a shape dependent friction along the whole trajectory, as it was done in HICOL, 
reproduces in a self-consistent way the mass distributions (Fig|](a)), the mass-angle corre- 
lations (Figs. H,0) and the fission fragment TKE. The other approaches discussed here, 
fail to reproduce all observables at once. There are minor discrepancies between the model 
and the experimental results which show that the form factor mixing window-plus-wall and 
mono-wall friction requires further investigation. The mass variances for the deep inelas- 
tic trajectories also need further refinements in order to yield better agreement with the 
experimentally measured cross-sections. 

E. Mass-Angle Correlations and Reaction Times 

From the measured angular distributions extrapolated to 0° and 180°, the double differ- 
ential cross-sections d 2 a / '(d9 c . m dA) were constructed. 
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The results are plotted in Fig. [11]. The contour levels are listed in the figure. A forward- 
backward asymmetry, which reveals the short time scale of the reaction, is clearly evident for 
both energies measured. A two-body friction mechanism with 1/T 2 temperature dependence, 
as suggested in Ref. [O], would require friction scaling factors of k° ut = 10 — 12 for E/ a b = 



383 MeV and k° ut = 4 for the higher energy, which in both cases results in time scale 
much longer than the rotational period of the system and flat da/d8 c , m , distributions. Only 
the deep inelastic component of the touch cross-section would retain the asymmetry about 
Oc.m. = 90°. Since the data do not show this behavior, the large friction scenario is ruled out. 
A direct comparison of the measured double differential cross-sections d 2 a j{dQ c , m dA) and 
HICOL could not be done, since the code predicts the mean value of the scattering angle, 
but does not calculate its variance. 

For all partial waves, the calculated mean scattering angles are smaller than 180°, which 
is consistent with the experimental observation of a non-orbiting angular distribution. Al- 
though HICOL suggests the correct behavior, there is an indication that the reaction times 
for the intermediate masses are over-predicted, since in the experiment, all fragments in this 
range are peaked forward with a maximum yield angle in the vicinity of 20° and an average 
angle of 60° — 80°. A larger saturation angle (« 100°) is obtained from the calculation. 
The correct description of the measured angular distributions would require a considerable 
dispersion around the mean values, which is not calculated by HICOL, and mean angles 
smaller by ~ 30°. 

An upper limit to the reaction time can be estimated using the characteristic rotational 
time of the system for the longest lasting trajectories. At E; a f, = 418 MeV, for the average 
quasi-fission trajectory, contributing in the symmetric mass bin A=120 — 140, L=46fr and 
moment of inertia of two touching spheres we estimate r ~ 30 x 10~ 21 s. In HICOL, the 
rotational frequency and the moment of inertia are changing along the trajectory and reflect 
the evolution of the shape and rotational degrees of freedom. HICOL's result for the same 
trajectory is r = 15 x 10~ 21 s, which is in excellent agreement with the time extracted from 
the neutron measurement, r = 18 ± 9 x 10~ 21 s 0. In contrast Wilczynska's calculation 
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13, extracted r = 80tf 5 x 10 21 s from the same data. Table IXTTI gives the HICOL results 



obtained for the average trajectories contributing in symmetric quasi-fission for the two 
energies measured. 

The two completely different physics approaches of Ref. || and this work have yielded 
similar results for the reaction time scale, which indicates that the one-body friction imple- 
mented in HICOL gives a satisfactory description of the reaction mechanism. However, a 
dynamical calculation of the pre-scission particle emission is necessary, in order to relate the 
pre-scission neutron multiplicity to the viscosity. Such calculations were done in this work 
and are presented in Appendix A. After adjusting the statistical model parameters, so that 
both the measured pre-scission neutron multiplicities and mean neutron kinetic energies are 
reproduced, the time scale derived from our dynamical calculation is in agreement with the 
original "static" results of Ref. |J. 

V. ONE-BODY DISSIPATION LIMITS 

To determine the experimental limits to the one-body viscosity, implemented in HICOL, 
calculations were done using scaling factors to the friction tensor ranging from k s = 0.5 to 
k s = 4 and tested versus the experimental mass and angular distributions. These factors 
were kept constant along the whole trajectory, thus preserving the original shape dependence 
of friction, which would be destroyed if different scaling factors were used in different parts 
of the trajectory. 



The results of the calculations are plotted in Fig.[12|. Scaling factors k s > 1.7 lead to 
scattering angles > 180° and orbiting type angular distributions, which are inconsistent with 
the data. In the lower limit, k s < 0.8, there is considerable discrepancy between the calcu- 
lated and the measured mass distributions. This analysis sets the upper and lower limits to 
the one-body friction scaling factors for the temperature range of the quasi-fission reactions 
studied in this work. The two measurements in the present study cover the temperature 
range of the pre-scission neutron data from Ref. M. In contrast to the calculations in Ref. 
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11] which found that scaling factors k s = 4 — 12, varying with temperature, were needed 



to describe the data in Ref. [[J, we consider both our measurements are consistent with 
one-body viscosity. Supporting this conclusion, Feldmeier in J7J analyzed the 238 U+ 48 Ca at 
Ezat> =5.9 MeV/u data from Shen et.al. ||, and reproduced the measured mass distributions 
with his code. A pre-scission GDR 7-ray measurement and analysis |25| in the quasi-fission 
reaction of 58 Ni+ 165 Ho at E/ a b = 368 MeV was also consistent with one-body friction. 

There have been a number of studies in fusion-fission and fusion-evaporation reactions 
based on the statistical model, which show the need for dissipation in the fission degree of 



freedom. The systematic studies by the Stony Brook group (See Ref. |T2[ and references 
therein) on pre-scission GDR 7— rays indicates that viscosity increases with temperature and 



is larger than the one-body value. Some of these results were used in Ref. [|ll[] and together 
with the pre-scission neutron multiplicity calculations discussed above were interpreted as 
evidence for strong two-body viscosity setting in at temperature T < 2 MeV. 

On the other hand, other studies from pre-scission particle multiplicities (e.g. Ref. [p6| , p7[ ) 
and GDR 7-rays(e.g.Ref. |28[]) in the same temperature range do not find evidence for such 
large viscosity. A recent measurement of Morton et.al. of GDR 7-ray multiplicities for 
the 32 S+ nai W, 208 Pb reactions also found a delay time consistent with one-body dissipation 
and 7-fission angular correlations that could be described without temperature dependent 
viscosity. At present, the dissipation mechanisms in fusion-fission reactions and the temper- 
ature at which the two-body dissipation sets in is still a subject of considerable debate. 



VI. CONCLUSIONS 

The present study concentrated on the quasi-fission reaction dynamics, exploiting the 
unique observables available in this type of reactions: mass-angle correlations and mass 
distributions. It aimed to clarify and reconcile the data and the theoretical descriptions of 
nuclear viscosity that were available prior to this. We find that two different experimental 
observables: mass-angle correlations, pre-scission neutron multiplicities have yielded similar 
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results for the reaction time scales. Our data analysis in terms of the one-body dissipation 
theory was qualitatively successful. Our own dynamical evaporation calculation with the 
code DIFHEAT was performed for the neutron data (See Appendix A) and reconciled with 
the original statistical model analysis of Hinde et.al.. In this work it was shown, how the 
mass-distributions by themselves are also sensitive to the dissipation mechanism. The new 
results do not confirm the onset of two-body viscosity in the temperature range below T=2.5 
MeV. 
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APPENDIX A: HICOL PLUS EVAPORATION 



In the Wilczynska et.al. [KJ approach, HICOL's trajectory was divided into time steps. 
A statistical model calculation was applied for each step, replacing the excitation energy 
with HICOL's dissipated energy minus the energy lost in evaporation from the previous 
step. 

A similar code, was developed independently in this work using a somewhat different 
approach. Instead of dividing the trajectory into time steps, a continuous evaporation 
process was assumed, following Refs. pO|j3T[] . The results of Refs. P!D]JII| could be reproduced 



with our code DIFHEAT. 

Calculating the accumulated neutron multiplicity along the L=40h trajectory in 
64 Ni+ 197 Au at Ei a b =418 MeV, using a level density parameter a = A/9 MeV -1 (the value 
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used in the neutron measurement), one finds that the calculated pre-scission multiplicity does 



not reach the experimental value. In Refs. |p!0[P!T| , the authors concluded that the reaction 
time predicted by HICOL using one-body friction is too short to be able to accommodate 
the measured pre-scission neutron multiplicities. 

A different approach to reconcile the measured pre-scission neutron multiplicity with the 
dynamical trajectory calculation was applied in this work. It appears that the discrepancy 



between Hinde's et.al. || and Wilczynski's et.al. jlQUTll calculations does not stem from the 
"static" versus "dynamical" approach, but rather reflects the different nuclear temperature 
achieved in the two cases. In Hinde's calculation the level density parameter was kept fixed 
at the value a = A/9 MeV -1 , but the excitation energy of the system was treated as a 
free parameter[| and adjusted in order to reproduce both the measured pre-scission neutron 
multiplicity and spectral shapes. Although routinely in statistical model calculations the 
level density parameter is varied, while the excitation energy is fixed, as Hinde pointed out, 



both observables are sensitive to the nuclear temperature T = ^jE*/a, which contains the 
ratio of the two quantities. The authors of Ref. [|Tt]] argued that the dynamically calculated 
excitation energy, as predicted by HICOL [] should be used. Keeping the level density 
parameter at the same value of a = A/9 MeV" 1 , but using lower excitation energy along 
the trajectory (E* ~ 100 MeV) results in lower nuclear temperature. Consequently, this 
calculation yields longer neutron lifetimes and, in addition, does not reproduce the mean 
neutron kinetic energy. 



The dashed curves in Fig. [13] show the neutron multiplicity and the mean neutron kinetic 
energy along the trajectory, calculated with a = A/ 9 MeV -1 and an input heat curve 
obtained by scaling the friction in HICOL by a factor k s = 4. The shaded regions indicate 



1 E* =E X + AEj;. ~E X = 95.3 MeV was determined from the Q value of the reaction and A E x = 50 
MeV was obtained from the fit to the data. 
2 See the total heat curve in Fig.||(c). 
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the experimental multiplicity and mean kinetic energy values reported in Ref. 0. 

Since HICOL indeed provides a more consistent description of the excitation energy, 
than the fitting procedure used in Ref. ||, in order to produce the same temperature (T« 
2.2 MeV), as derived from the neutron spectra measured by Hinde, we need to vary the 
level density parameter. From Fig. ^(c) we find that E* ~ 100 MeV, and we estimate that 
a = A/12.7 MeV -1 should reproduce both the neutron lifetimes and mean kinetic energy. 
The solid lines in the Fig. [13] show the result from DIFHEAT, using a level density parameter 
a = A/13 MeV -1 and a one-body dissipation heat curve coming from HICOL. 

Unlike the conclusion in Ref. ||11|| , that large friction is necessary to reproduce the exper- 
imental pre-scission neutron multiplicities, we find that the reaction time extracted from the 
neutron data ~ 20 x 10~ 21 s is consistent with HICOL's predictions and with the evaporation 
calculation using DIFHEAT, if the statistical model parameters are properly constrained. 
While HICOL is free of adjustable parameters, in any statistical model calculation one needs 
at least two observables, to be able to constrain the multi-parameter model calculation. It is 
beyond the scope of this work to investigate the influence of all statistical model parameters 
on the results from DIFHEAT. There is a significant deformation dependence both in the 
particle transmission coefficients and their binding energies P2"|j2"E|| . Accounting for these 
effects might change the value of the level density parameter needed to describe the data. 
However, we note that values as small as a ~ A/13 MeV -1 have been previously deduced 
from particle spectral shapes in heavy hot systems |]33| , |34 . 
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TABLES 



Reaction 



E lab MeV 



Tf X 10~ 21 S a TfX W~ 21 S b 



40 Ar+ 23 8u 



249 



3.25±0.2 



30±10 



390± 



60 



64 M+ 197 Au 



418 



3.15±0.6 



18±9 



80±H 



64 Ni+ 208p b 



418 



3.25±0.6 



28±10 



70±|| 



64 Ni+ 238u 



418 



4.00±0.8 



15±8 



120±18 



TABLE I. Quasi-fission reactions, pre-scission neutron multiplicities and reaction times. 
a from Ref. g 



Trom Ref. p| 
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E lab = 418 MeV 


E lab = 383 MeV 


Mass 
( arnn I 


cr c 
(mb) 


(mb) 


@touch 

(mb) 


(mb) 


(mb) 


&touch 

(mb) 


120-180 


95 ±21 




95 ±21 


43 ± 10 




43 ± 10 


110-120 


128 ± 29 




128 ± 29 


46 ± 11 




46 ± 11 


100-110 


162 ± 37 




162 ± 37 


63 ± 15 




63 ± 15 


90-100 


199 ± 45 




199 ± 45 


70 ± 16 




70 ± 16 


70-80 
60-70 
50-60 


onfi 4- A 7 

1 37+ 34 
J-" 3 ' -68 

9512! 
291L 


7 c+38 
' °-19 

186li 
11311 


onfi 4_ A 7 
ZUD ± 4 I 

212 ± 59 
281 ± 79 
142 ± 40 


1 net 4_ O/i 
1UD ± Z4t 

501" 
6.3l 3 i 


62t£ 

1941S 

134+ 6 33 7 


i net 4- 9/1 

1UD ± Z4t 

124 ± 32 
245 ± 65 
140 ± 36 


total 




374+$° 


1427 ± 357 






841 ± 210 



TABLE II. Angle-integrated cross-sections in bins of 10 mass units from the two energies 
studied. 
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E lab (MeV) 


(L)h 


time x 10 21 s 


418 


46 


15 


383 


22 


16 


TABLE III. 


Reaction times for the average 


quasi-fission trajectories contributing in the 



A= 120 — 140 mass range. 
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FIG. 1. The experimental layout in the "BigMac" scattering chamber. 
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FIG. 2. Scatter plots from detector 3 located at 40° in the lab system for E; a b = 418 MeV. (a) 
Energy versus TOF (uncalibrated). The "banana" gate applied to select the events of interest is 
shown, (b) TKE (MeV) versus Mass (amu). The gate shown separates the quasi-fission and the 
deep inelastic events from the elastic and quasi-elastic scattering. The line shows Viola systematics 
TKE. 
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FIG. 3. Efficiency curves for all 8 detectors in the two positions of the platform for E; a 5 = 418 
MeV and E lab = 383 MeV. 
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FIG. 4. Fragment angular distributions in bins of 10 mass units from the E; a b =418 MeV 
measurement. The fits to the data are described in the text. The values of the measured cross 
sections are given. The open points come from the measurements in the complementary mass bin 
with the assumption that binary fission is observed. 
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FIG. 5. As per Fig.| but for E lab = 383 MeV. 
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FIG. 6. Along the trajectory calculations: (a) potential energy and shapes. (b)Kinetic energy in 
the shape degrees of freedom T s , relative rotation T re i and intrinsic rotation Tj nt . (c) Dissipated 
energy: total heat and the relative contributions from the window and wall dissipation. 
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FIG. 7. Window-plus-wall or mono-wall? (a) The mass evolution of the projectile-like and the 
target-like fragments is calculated for the Ei a b = 418 MeV, L=50?i trajectory with window-plus-wall 
and mono- wall friction mixed with a from factor f{x n ) as described in the text, (b) The evolution 
of the form factor f(x n ) along the same trajectory, (c) Fragment mass evolution calculated without 
the window friction ( f(x n ) = 1 along the whole trajectory independent of the shape) and (d) the 
mass evolution with f(x n ) = - window-plus-wall everywhere. 
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FIG. 8. As per Fig.^ but calculated with large friction as in Ref. [11]. 
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FIG. 9. Mass distributions: (a) E ia6 = 418 MeV and (b) Ei ah = 383 MeV. The data are shown 
with points. The histograms give HICOL calculations: dashed lines - unmodified HICOL; solid 
lines - the window-plus-wall to mono-wall form factor is modified as described in the text. 
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FIG. 10. Mass distribution calculations for = 418 MeV in comparison to the data: (a) 
a calculation without the window friction; (b) a calculation with friction scaling according to 
Wilczynski; (c) 0.3 x (window+wall) along the whole trajectory. 
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E=418 MeV , Mass-Angle Correlations^ 
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FIG. 11. Double differential cross-sections d 2 a / {d9 cm dA) (mb/rad/amu) from the two energies 
measured. The contour levels are shown in the plots. 
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FIG. 12. Mass distributions (a) and mean scattering angles (b) calculated with HICOL with 
scaling factors : ks = 1.7 - long-dashed histograms and curves, ks = 0.8 - short-dashed histograms 
and curves, ks = 1 - full lines. The experimental data for = 418 MeV are shown with points. 



40 



5 
4 
3 
2 





5 

4 
3 
2 







i i r 



i i i i i r 



n r 



64 Ni + i97 Au E lab =418 MeV, L=40h 



experimental multiplicity 




k =1 



k =4 



a=A/13 MeV-i ^ a=A/9 MeV -1 



= i i r 



"i i i i i r 



experimental {^ kiD ) 







200 400 600 

time x 10 _22 s 



i - 



i i i i i i i i i i i i i i i 



800 



13. Pre-scission neutron multiplicity and mean neutron kinetic energy along the trajec- 
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